BELLE-CONF-0743 



Measurement of B —> D s Kn branching ratios 

K. Abe, 10 I. Adachi, 10 H. Aihara, 52 K. Arinstein, 1 T. Aso, 56 V. Aulchenko, 1 T. Aushev, 22 - 16 
T. Aziz, 48 S. Bahinipati, 3 A. M. Bakich, 47 V. Balagura, 16 Y. Ban, 38 S. Banerjee, 48 
E. Barberio, 25 A. Bay, 22 I. Bedny, 1 K. Belous, 15 V. Bhardwaj, 37 U. Bitenc, 17 S. Blyth, 29 
A. Bondar, 1 A. Bozek, 31 M. Bracko, 24 ' 17 J. Brodzicka, 10 T. E. Browder, 9 M.-C. Chang, 4 

P. Chang, 30 Y. Chao, 30 A. Chen, 28 K.-F. Chen, 30 W. T. Chen, 28 B. G. Cheon, 8 
C.-C. Chiang, 30 R. Chistov, 16 I.-S. Cho, 58 S.-K. Choi, 7 Y. Choi, 46 Y. K. Choi, 46 S. Cole, 47 
J. Dalseno, 25 M. Danilov, 16 A. Das, 48 M. Dash, 57 J. Dragic, 10 A. Drutskoy, 3 S. Eidelman, 1 
D. Epifanov, 1 S. Fratina, 17 H. Fujii, 10 M. Fujikawa, 27 N. Gabyshev, 1 A. Garmash, 40 
A. Go, 28 G. Gokhroo, 48 P. Goldenzweig, 3 B. Golob, 23 ' 17 M. Grosse Perdekamp, 12 ' 41 
H. Guler, 9 H. Ha, 19 J. Haba, 10 K. Hara, 26 T. Hara, 36 Y. Hasegawa, 45 N. C. Hastings, 52 
K. Hayasaka, 26 H. Hayashii, 27 M. Hazumi, 10 D. Heffernan, 36 T. Higuchi, 10 L. Hinz, 22 
H. Hoedlmoser, 9 T. Hokuue, 26 Y. Horii, 51 Y. Hoshi, 50 K. Hoshina, 55 S. Hou, 28 
W.-S. Hou, 30 Y. B. Hsiung, 30 H. J. Hyun, 21 Y. Igarashi, 10 T. Iijima, 26 K. Ikado, 26 
K. Inami, 26 A. Ishikawa, 42 H. Ishino, 53 R. Itoh, 10 M. Iwabuchi, 6 M. Iwasaki, 52 Y. Iwasaki, 10 
C. Jacoby, 22 N. J. Joshi, 48 M. Kaga, 26 D. H. Kah, 21 H. Kaji, 26 S. Kajiwara, 36 
H. Kakuno, 52 J. H. Kang, 58 P. Kapusta, 31 S. U. Kataoka, 27 N. Katayama, 10 H. Kawai, 2 
" Kawasaki, 33 A. Kibayashi, 10 H. Kichimi, 10 H. J. Kim, 21 H. 0. Kim, 46 J. H. Kim/ fi 
K. Kim, 44 Y. J. Kim, 6 K. Kinoshita, 3 S. Korpar, 24 ' 17 Y. Kozakai, 26 P. Krizan, 23 * 
. Krokovny, 10 R. Kumar, 37 E. Kurihara, 2 A. Kusaka, 52 A. Kuzmin, 1 Y.-J. Kwon, f 

C T 5 n T „A™ 14 T T„„44 7 c T 46 T T 44 Q T? T ™ 44 rp T • i ; 



r. iyi uisjjv iiy , rt. ivuiiicii, hi. rvuniiciici, i\. ivuscu^ci, i\. ivuzimii, i.-j. ivwuu, 

J. S. Lange, 5 G. Leder, 14 J. Lee, 44 J. S. Lee, 46 M. J. Lee, 44 S. E. Lee, 44 T. Lesiak, 31 
J. Li, 9 A. Limosani, 25 S.-W. Lin, 30 Y. Liu, 6 D. Liventsev, 16 J. MacNaughton, 10 
G. Majumder, 48 F. Mandl, 14 D. Marlow, 40 T. Matsumura, 26 A. Matyja, 31 S. McOnie, 47 
T. Medvedeva, 16 Y. Mikami, 51 W. MitarofT, 14 K. Miyabayashi, 27 H. Miyake, 36 H. Miyata, 33 

Y. Miyazaki, 26 R. Mizuk, 16 G. R. Moloney, 25 T. Mori, 26 J. Mueller, 39 A. Murakami, 42 
T. Nagamine, 51 Y. Nagasaka, 11 Y. Nakahama, 52 I. Nakamura, 10 E. Nakano, 35 M. Nakao, 10 
H. Nakayama, 52 H. Nakazawa, 28 Z. Natkaniec, 31 K. Neichi, 50 S. Nishida, 10 K. Nishimura, 9 
Y. Nishio, 26 I. Nishizawa, 54 0. Nitoh, 55 S. Noguchi, 27 T. Nozaki, 10 A. Ogawa, 41 
S. Ogawa, 49 T. Ohshima, 26 S. Okuno, 18 S. L. Olsen, 9 S. Ono, 53 W. Ostrowicz, 31 
H. Ozaki, 10 P. Pakhlov, 16 G. Pakhlova, 16 H. Palka, 31 C. W. Park, 46 H. Park, 21 
K. S. Park, 46 N. Parslow, 47 L. S. Peak, 47 M. Pernicka, 14 R. Pestotnik, 17 M. Peters, 9 
L. E. Piilonen, 57 A. Poluektov, 1 J. Rorie, 9 M. Rozanska, 31 H. Sahoo, 9 Y. Sakai, 10 
H. Sakamoto, 20 H. Sakaue, 35 T. R. Sarangi, 6 N. Satoyama, 45 K. Sayeed, 3 T. Schietinger, 22 
O. Schneider, 22 P. Schonmeier, 51 J. Schumann, 10 C. Schwanda, 14 A. J. Schwartz, 3 
R. Seidl, 12 ' 41 A. Sekiya, 27 K. Senyo, 26 M. E. Sevior, 25 L. Shang, 13 M. Shapkin, 15 
C. P. Shen, 13 H. Shibuya, 49 S. Shinomiya, 36 J.-G. Shiu, 30 B. Shwartz, 1 J. B. Singh, 37 
A. Sokolov, 15 E. Solovieva, 16 A. Somov, 3 S. Stanic, 34 M. Staric, 17 J. Stypula, 31 
A. Sugiyama, 42 K. Sumisawa, 10 T. Sumiyoshi, 54 S. Suzuki, 42 S. Y. Suzuki, 10 O. Tajima, 10 

F. Takasaki, 10 K. Tamai, 10 N. Tamura, 33 M. Tanaka, 10 N. Taniguchi, 20 G. N. Taylor, 25 
Y. Teramoto, 35 I. Tikhomirov, 16 K. Trabelsi, 10 Y. F. Tse, 25 T. Tsuboyama, 10 K. Uchida, 9 
Y. Uchida, 6 S. Uehara, 10 K. Ueno, 30 T. Uglov, 16 Y. Unno, 8 S. Uno, 10 P. Urquijo, 25 



Typeset by REVTeX 



1 



Y. Ushiroda, 10 Y. Usov, 1 G. Varner, 9 K. E. Varvell, 47 K. Vervink, 22 S. Villa, 22 
A. Vinokurova, 1 C. C. Wang, 30 C. H. Wang, 29 J. Wang, 38 M.-Z. Wang, 30 P. Wang, 13 
X. L. Wang, 13 M. Watanabe, 33 Y. Watanabe, 18 R. Wedd, 25 J. Wicht, 22 L. Widhalm, 14 
J. Wiechczynski, 31 E. Won, 19 B. D. Yabsley, 47 A. Yamaguchi, 51 H. Yamamoto, 51 
M. Yamaoka, 26 Y. Yamashita, 32 M. Yamauchi, 10 C. Z. Yuan, 13 Y. Yusa, 57 C. C. Zhang, 13 
L. M. Zhang, 43 Z. P. Zhang, 43 V. Zhilich, 1 V. Zhulanov, 1 A. Zupanc, 17 and N. Zwahlen 22 

(The Belle Collaboration) 

1 Budker Institute of Nuclear Physics, Novosibirsk 

2 Chiba University, Chiba 
3 University of Cincinnati, Cincinnati, Ohio 45221 
4 Department of Physics, Fu Jen Catholic University, Taipei 
5 Justus-Liebig-Universitat Giefien, Giefien 
6 The Graduate University for Advanced Studies, Hayama 
7 Gyeongsang National University, Chinju 
8 Hanyang University, Seoul 
9 University of Hawaii, Honolulu, Hawaii 96822 
10 High Energy Accelerator Research Organization (KEK), Tsukuba 
11 Hiroshima Institute of Technology, Hiroshima 
12 University of Illinois at Urbana- Champaign, Urbana, Illinois 61801 
13 Institute of High Energy Physics, 
Chinese Academy of Sciences, Beijing 
^Institute of High Energy Physics, Vienna 
15 Institute of High Energy Physics, Protvino 
16 Institute for Theoretical and Experimental Physics, Moscow 
17 J. Stefan Institute, Ljubljana 
18 Kanagawa University, Yokohama 
19 Korea University, Seoul 
20 Kyoto University, Kyoto 
21 Kyungpook National University, Taegu 
22 Ecole Polytecnique Federale Lausanne, EPFL, Lausanne 
23 University of Ljubljana, Ljubljana 
2J " University of Maribor, Maribor 
25 University of Melbourne, School of Physics, Victoria 3010 
26 Nagoya University, Nagoya 
27 Nara Women's University, Nara 
28 National Central University, Chung-li 
29 National United University, Miao Li 
30 Department of Physics, National Taiwan University, Taipei 
31 H. Niewodniczanski Institute of Nuclear Physics, Krakow 
32 Nippon Dental University, Niigata 
33 Niigata University, Niigata 
3J ' University of Nova Gorica, Nova Gorica 
35 0saka City University, Osaka 

36 0saka University, Osaka 
37 Panjab University, Chandigarh 
38 Peking University, Beijing 



2 



39 University of Pittsburgh, Pittsburgh, Pennsylvania 15260 

4° Princeton University, Princeton, New Jersey 08544 
41 RIKEN BNL Research Center, Upton, New York 11973 
42 Saga University, Saga 
^ 3 University of Science and Technology of China, Hefei 
^ Seoul National University, Seoul 

^ 5 Shinshu University, Nagano 
^ Sungkyunkwan University, Suwon 
47 University of Sydney, Sydney, New South Wales 
^ 8 Tata Institute of Fundamental Research, Mumbai 
^ 9 Toho University, Funabashi 
50 Tohoku Gakuin University, Tagajo 
51 Tohoku University, Sendai 
5i ''Department of Physics, University of Tokyo, Tokyo 

53 Tokyo Institute of Technology, Tokyo 

54 Tokyo Metropolitan University, Tokyo 

55 Tokyo University of Agriculture and Technology, Tokyo 
56 Toyama National College of Maritime Technology, Toyama 
57 Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061 

58 Yonsei University, Seoul 

Abstract 

We report a measurement of the exclusive B + meson decay to the final state D~K + tt + using 
520 x W 6 BB pairs collected near the T(4S) resonance, with the Belle detector at the KEKB 
asymmetric-energy e + e~ collider. Using the Dg — > 4>ir~ decay mode to reconstruct D~ mesons, we 
obtain the branching fraction B{B + -► D~K + ir + ) = (1.77±%\l(stat)±0.16(syst)±0.23(B)) x 1(T 4 . 
We also present preliminary results of a study of the two-body D S K, D s ir and Kit subsystems 
observed in this final state. 

PACS numbers: 13.25.Hw, 14.40.Nd 
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I. INTRODUCTION 



We search for the exclusive decays of charged B mesons into D s Kn final states. These 
modes offer rich possibilities for studies of different two-body subsystems, such as D S K. For 
the first decay mode studied, B + — > D~ K + n + [1], the dominant process is described by the 
Feynman diagram shown in Fig. 1. This process is mediated by the b —>■ c quark transition 
and includes the production of an ss pair via a radiative gluon. The D s Ktt final state can 
also be a result of the b — > c — > s decay chain, where the D s meson and the charged pions are 
produced in two W boson decays. The Feynman diagram in Fig. 2 describes the dominant 
process responsible for the second decay mode, B + — > D+D° with D° — > K + tt~ . Note that 
the two different processes lead to a similar three-body final state, but with opposite charges 
for the D s and 7r mesons. We measure the branching fraction of the first decay mode and 
study two-body subsystems, while the second decay mode is a control channel to check the 
reliability of the method. 

In the following the Belle detector and the data sample are briefly described. Next 
we discuss the reconstruction of the D s Ktt final states. Finally the determination of the 
branching fractions together with a preliminary study of the D s Kn subsystems including 
D S K, D S 7T and Kit are presented. 




FIG. 1: Diagram for the decay B + — * D s K + ir 




FIG. 2: Diagram for the decay B+ -» D+D°, D° -> K+ir- . 
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II. DETECTOR AND DATA SAMPLE 



The results are based on a data sample that contains (520.2 ± 6.8) x 10 6 BB pairs, cor- 
responding to an integrated luminosity of 479 fb -1 , collected with the Belle detector at the 
KEKB asymmetric-energy e + e~ (3.5 on 8 GeV) collider [2]. KEKB operates at the T(4S f ) 
resonance = 10.58 GeV) with a peak luminosity that exceeds 1.7 x 10 34 cm~ 2 s _1 . The 
T(4S) resonance is produced with a Lorentz boost of (5^ = 0.425 nearly along the electron 
beamline (z-axis). The production rates of B + B~ and B°B° pairs are assumed to be equal. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter composed of CsI(Tl) crystals (ECL) located 
inside a super-conducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect K\ mesons and to identify muons 
(KLM). The detector is described in detail elsewhere [3]. Two inner detector configurations 
were used. A 2.0 cm beampipe and a 3-layer silicon vertex detector was used for the first 
sample of 152 x 10 6 BB pairs, while a 1.5 cm beampipe, a 4-layer silicon detector and a 
small-cell inner drift chamber were used to record the remaining 368 x 10 6 -B-B pairs [4]. 

III. RECONSTRUCTION OF B+ -► DjK + ir + AND B+ -► D+D° 

The reconstruction of B + — > D~K + n + and B + — > DfD° decays consists of the following 
steps: selection of charged tracks, discrimination between BB and continuum events, particle 
identification and reconstruction of all intermediate decays in B — > D s Kn processes. Each 
of these steps is briefly described below. 

A. Track selection 

Charged tracks with dr < 0.5 cm and dz < 5 cm are selected, where dr and dz are impact 
parameters measured in the r-<f> (transverse) plane and z direction, respectively. Charged 
tracks are also required to have transverse momenta greater than 100 MeV/c. 

B. Suppression of continuum events 

We exploit the event topology to discriminate between spherical BB events and the dom- 
inant background from jet-like continuum events, e + e~ — > qq (q = u, d, s, c). We combine 
event shape variables using all particles in an event to calculate Fox- Wolfram moments [5] 
and define the R 2 , ratio of second and zeroth moments, as 

= Eij \Pi\\Pj\P2(cO8 0ij) 

T,i,j\Pi\\Pj\ 

where p indicates the particle momentum, P 2 is the Legendre polynomial of second order, 
and i, j enumerate all particles in the event. We require R 2 to be less than 0.4. 
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C. Particle identification 



Hadron identification is based on information provided by the CDC, ACC and TOF. For 
kaons we require Ck/it > 0.6 and C p /k < 0.95 (veto), where C x / y = c c ^ c (x, y — n, K, p) 
denotes the corresponding likelihood ratio. Pions are selected as non-kaons, satisfying veto 
conditions for K and p: C p /k < 0.95 and Cr/w < 0.95. In addition, we reject tracks that are 
consistent with an electron hypothesis. A selection imposed on this ratio results in a typical 
kaon (pion) identification efficiency ranging from 92% to 97% (94% to 98%) for various 
decay modes, while 2% to 15% (4% to 8%) of kaon (pion) candidates are misidentified pions 
(kaons) . 

D. Reconstruction of exclusive Df decays 

The candidates are reconstructed in two final states: 0(— > K + K~)ii + and 
i?*(892)°(^ K-it + )K + . We reconstruct <\> mesons in the K + K final state. The 6 mass 
is estimated to be (1019.70 ± 0.153) MeV/c 2 , in agreement with the world average value 
mwA{4>) = (1019.460 ± 0.019) MeV/c 2 [6]. We accept K + K~ pairs satisfying the require- 
ment \m(K + K~) — mwA{4>)\ < 10 MeV/c 2 . Similarly, the K + n~ mass spectrum exhibits a 
7^(892)° signal at a mass of (892.19 ± 1.80) MeV/c 2 {m WA (K* Q ) = (896.00 ± 0.25) MeV/c 2 
[6]); K*° candidates are thus selected with the requirement: \m(K + ir~) — m WA (K*°)\ < 100 
MeV/c 2 . 

A clear D s meson signal is observed in both channels considered (Figs. 3-4). The 
fitting procedure - described in the next section - is used to obtain the mean values of 
the D s mass in each channel. This mass, averaged over two decay channels studied, is 
estimated to be (1968.1 ± 0.20(stat)) MeV/c 2 , which is consistent with the world average 
value i71wa{D s ) = (1968.2 ± 0.5) MeV/c 2 [6]. The invariant mass of D s candidates is 
required to satisfy the criterion \m(D s ) — my/ a{D s )\ < 15 MeV/c 2 , corresponding to a 
window of approximately three standard deviations about the D s mass. 

E. Reconstruction of B mesons 

B mesons are reconstructed combining D s candidates with identified kaons and pions. In 
exclusive reconstruction of B mesons two kinematic variables are used: the energy difference, 
AE, and the beam-energy-constrained mass, M bc . These are defined as: 

AE = E B — -Ebeam (1) 

and 

M^ = yjEL m -p l B , (2) 

where E B and P B are the reconstructed energy and momentum of the B candidate, and 
-Kbeam is the run-dependent beam energy, all expressed in the centre-of-mass (CM) frame. 
We use the dedicated Monte Carlo (MC) sample of B — > D s Kn decays to define signal 
regions in AE and M bc as \AE\ < 0.03 GeV and M bc > 5.27 GeV/c 2 . For all B candidates, 
a loose requirement on the goodness of the D s Kn vertex fit - with the D s mass constrained 
to the world average value - is also applied (x 2 b < 50). 
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From MC simulation we determine that the background contribution from B + — > ((cc) — > 
i^ + i^ _ 7r + 7r _ )i^ + decays, where (cc) are charmonium states such as the J/vp or r] c , is removed 
by discarding all events that satisfy the requirement: 2.88 GeV/c 2 < M (K + K~7i + 7i~) < 
3.18 GeV/c 2 . 

For B + — > DfD° candidate events a requirement on fT + 7r~ invariant mass is imposed: 
\m(K + n^)— m W A(D°)\ < 0.015 GeV/c 2 [6]. This condition is determined by using DfK + n~ 
combinations for events from the (AE, M^ c ) signal region and plotting the invariant mass 
of K + ir~ pairs as shown in Fig. 5. 

After selection requirements, 8.8%(10.6%) of events with reconstructed D s — > (f)ir (D s — > 
K*°K) decays have more than one B candidate. In these cases we select the B candidate 
having the smallest value of \b- If there are two or more combinations with the same value 
of Xb, the one containing a kaon - originating directly from the B decay - with the highest 
likelihood ratio C K / W is selected. 

For events passing all selection criteria mentioned above the AE distributions (requiring 
Mb c > 5.27 GeV/c 2 ) and M^ c distributions (requiring \AE\ < 0.03 GeV) are examined. The 
plots shown in Figs. 6-9 exhibit clear signals from the corresponding B meson decays. In 
addition, an excess of events can be seen in the range AE e (—0.2, —0.08) GeV, probably 
due to events with one unreconstructed particle, such as a n° or 7, in the final state. 



IV. DETERMINATION OF B — > D s Kir BRANCHING FRACTION 



The branching fractions are calculated using the formula: 

B(B - D s Kn) = Ns,£ (3) 

ly BB ' °int 

where N BB - is the number of produced BB pairs, N$ is the number of reconstructed B — > 
D s Kn events, e denotes the reconstruction efficiency and Bi nt is the product of branching 
fractions of intermediate resonances, present in the respective B — > D s Kix processes (values 
summarised in Table I are taken from the Ref. [6]). 



TABLE I: Branching fractions of intermediate resonances present in B + — > D s K + tt + and B + — > 
D+D° decays (from the PDG). 



Decay 


Branching fraction [%] 




2.16 ± 0.28 


D+ -> K*(892)°K + , i?*(892)° -> K~tt + 


2.5 ± 0.5 




3.80 ± 0.07 



A. Determination of the number of B — > D s Kir events 



The yields of B — > D s Kn events are determined from an unbinned extended maximum 
likelihood fit to (AE, M bc , m Ds ) distributions using the MINUIT [7] package. The likelihood 
function is given by: 



N 



c(e) = ]Jf(x i 
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where x*i is a vector of independently measured values of the probability density function 
f(x\8). The index % {% — 1, 2, . . ., N) counts the number of reconstructed events used in the 
fit. The vector 6 = (9\, . . . ,9 n ) contains the parameters that are obtained by maximizing 
the likelihood function: 

9 = max£(#) 

S 

The probability density function / contains a Poisson factor to include the relation between 
the estimated number of all events (N s + N B ) and the number of reconstructed events N, 
and is defined as: 

f(x\8) = g(x\9) ■ ±-(N 8 + N B ) N e~ Ns ~ Ne , 



while 

g{x\9) 



N s 



N S + N B a AE V2n 



exp 



AE 2 



a M , V2n 



exp 



(M bc - m WA {B)f 



V2tt \ a, 



■ exp 



(m Ds - m WA (D s )Y 



2a, 



m_D„i 



+ 



■ exp 



{m Ds - m WA (D s )Y 
In 2 




iAE + w 2 AE 2 ) ■ 



f M hc y 
Kvm) exp 



yTs/2 

(po + pim Ds + p 2 m Ds 2 ) . 



arg 



(a) 
(b) 

) 

(c) 
(d) 

(e) 
(f) 



The first three contributions in this formula ((a),(b) and (c)) describe the signal parame- 
terized by four Gaussian functions: one for AE, one for Mb c and two Gaussians with the 
same mean and different widths for the rriD s distributions. The last three terms describe 
the background: second-order polynomials for AE and rri£, s variables ((d) and (f)), and the 
so-called ARGUS function [8] in (e) for the M^ c variable. 

In this notation x = (AE, M hc ,m Ds ) and 9 = (N s ,N B ,wo,Wi,W2,Po,Pi,P2,arg,<TAE, 
a Mhc ,a rnDsl ,(7 mDs2 ). Here N s and N B denote the number of signal events and the number 
of background events, respectively, while wo, wi,W2,po,Pi,P2 are the parameters of the poly- 
nomials, arg is a parameter of the ARGUS function and <jae, °"M bc , % sl , cr m Dg2 are the 
respective widths of the Gaussian functions. The factor M provides a normalization. 



TABLE II: Values of the Mbc AE and mo s parameters, as extracted from the fits described in 
the text. 



Decay 


M bc [MeV/c 2 ] 


AE [MeV] 


m Ds [MeV/c 2 ] 


B+ -► 




4>tt )K + t: + 


5279.4 ± 0.2 


-1.3 ±0.7 


1967.8 ± 0.3 


B+ -► 


DA^ 


K* (892)° K-)K + ir + 


5279.2 ±0.2 


-2.1 ±0.9 


1968.7 ±0.4 


B+ 


Df(^ 




5279.2 ±0.1 


-0.3 ±0.5 


1968.2 ± 0.2 


B+ -► 




K*(892)°K+)D° 


5279.1 ±0.1 


-0.2 ±0.6 


1967.9 ± 0.2 
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When performing the fit, all 9 parameters are allowed to float. The mean values of 
AE, Mbc and rriD s obtained from the fit are collected in Table II. The quoted M^ c (j/idJ 
values are mutually consistent and also in agreement with the world average value mwA(B) 
(m WA (D s )) [6]. 

The determination of the number of reconstructed decays (N) proceeds in two steps: 

• The fit to the (AE, M^ c , m,D s ) distribution is evaluated in the signal region (with 
requirements on all intermediate resonances), resulting in N events. 

• For decays containing K*(892)° meson, a fit to the sample of events from the if* (892)° 
mass sidebands ([0.746,0.796] GeV/c 2 and [0.996,1.046] GeV/c 2 ) is also performed. It 
yields the number of background events, N±, peaking in the signal region and therefore 
contributing to the iVo value. The same procedure is applied also for the efficiency 
calculation. This step is motivated by a MC study, where such a background contri- 
bution was found. The study also revealed that for the D s — * 0rr decay mode this 
background contribution is negligible. 

Finally, the number of B — > D s Kn events is calculated as: 

N = N - JVi. (4) 

The values of iV , Ni and N for the analysed decays are listed in Table III. 



TABLE III: Numbers of the signal and background events for analysed decays and efficiencies of 
their reconstructions. The errors are statistical only. 



Decay 



B + -► Ds((f)lT)K + TT + 

B+ -► D-(K*°K-)K + tt^ 
B+ -► D+((j)ir)D 
B+ -► D+(K*°K+)D° 



N 



N 



e 7o 



276.5 t\H 
299.5 ' 



512.3 



-26.0 
-25.1 
+23.4 
-22.8 



48.2 



+18.6 
-13.7 



276.5 ±™ 
251.3 ' 



512.3 



K32.0 
-28.6 
+23.4 
-22.8 



508.0 



47.8 



-10.4 
-7.6 



460.2 +Hi 



13.89 ±0.13 
9.05 ±0.11 
14.09 ±0.13 
10.03 ±0.12 



TABLE IV: Values of contributions (in %) to the overall systematic uncertainty in the branching 
fraction for the B + — > D~K + n + decay. 



Individual 


Values 


[%] for studied decays: 




contribution 


B+ 


-► D-((f)7r-)K + TT + 


B + D-(K*°K-)K+ir+ 


AB £ 




5.05 


4.16 


AB rn9 




1.26 


2.39 


AB r 




1.57 


8.04 


AB sel 




2.43 


0.72 


AB id 




5 


5 


^f^trck 




5 


5 


AB BB 




1.3 


1.3 


total 




9.34 


11.83 
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B. Determination of the reconstruction efficiency 

The reconstruction efficiency (e) of each channel in question is determined using MC 
samples of 10 5 e + e~ — > T(AS) — > B + B~ decays. These samples are generated in two subsets 
corresponding to the two inner detector configurations that were used in Belle experiment 
(denoted as SVD1 and SVD2). Each subset contains 5 x 10 4 events. In each event, one 
of the charged B mesons decays to the appropriate D : fK ± n ± (D+D°/D~D°) final state. 
These dedicated MC samples are subjected to the same analysis as the data, yielding the 
number of reconstructed events, NmCi and Nmc 2 i which corresponds to the SVD1 and 
SVD2 samples, respectively. These values divided by the number of all simulated events, 
gives partial reconstruction efficiencies for the studied decays: 

_ N M d _ N MC2 

£l ~ 50000' 62 ~ 50000 

The final efficiences for all studied decays are obtained by calculating the weighted averages 
of the values defined above: 

N BBl ■ e 1 + N BB2 ■ e 2 

€ = 

N BBl + N B B2 

where N BBl and N BB2 are the numbers of BB events collected during the SVD1 and SVD2 
running periods, respectively. Values of efficiencies for studied decays are presented in 
Table III. 

C. Studies of systematic uncertainties 

The following sources of systematic uncertainties are taken into account (cf. Table IV): 

• AB £ - uncertainty of efficiency determination estimated as a statistical error of its 
measurement, i.e. 

Ae = J( ^3*1 V . (A £l ) 2 + ( — V • (Ae 2 ) 2 

)j\N BBl + N BB2 ) { 11 \N BB1 + N BB2 ) 1 2) 

where Ae± and Ae 2 are statistical uncertainities on partial efficiencies. 
Considering the control channel and taking into account a mismatch between branch- 
ing fractions calculated for the MC sample (cf. IV. D subsection), an additional contri- 
bution to the efficiency systematics can be assigned. The shift (with an error) between 
the obtained branching fraction value and the generator value is calculated and the 
corresponding systematic error is determined. This factor is summed in quadrature 
with the above Ae value to evaluate the final AB e contribution. 

• AB sel - uncertainty due to the selection procedure. To estimate this error, the require- 
ment for the i?2 parameter is varied: R 2 < 0.40 — > R 2 < 0.35. 

• AB rng - error from changing the range of the fit to the AE variable. For B + — > 
D-K+U+ decays this range is changed from (-0.08,0.2) GeV to (-0.12,0.2) GeV. 
For each case the number of signal events iV™ s is determined, and the deviation from 
the iV value is calculated. 
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• AB T - uncertainty in the signal width. The width of the signal peak, as obtained from 
the fit to the AE/M^, C using MC sample, is then used as a fixed parameter in the same 
fit to data. The fixing procedure is applied sequentially: for AE, for M^ c and for both 
the AE and M bc signal peaks. Finally, the maximum deviation from the N value is 
chosen. 

• AB ld - uncertainty of particle identification in Belle experiment. A standard value 
of 1% per charged particle is assumed and uncertainties from kaons and pions are 
combined linearly, thus giving an overall contribution of 5 %. 

• A£>' rcfc - uncertainty of track reconstruction. A standard value of 1% per charged track 
is assumed giving an overall contribution of 5 %. 

• AB BB - uncertainty of the number of BB mesons used as data sample. Value of 1.3% 
is assumed in each studied decay. 

All contributions are assumed to be independent, hence the overall systematic error is ob- 
tained by summing those contributions in quadrature. The last three sources of systematic 
uncertainties are assumed to be common to all decays studied. All information about sys- 
tematic uncertainties is collected in Table IV. 



TABLE V: Branching fraction of studied decays. 



studied 


B 


uncertainties 


signif. 


decays 


xlO" 4 


stat.(+) 


stat.(-) 


syst. 


Bint 


M 


B + -► D-((f)7T)K + TT + 


1.77 


0.12 


0.12 


0.16 


0.23 


27.5 


B+ -► D-(K*°K-)K + tt + 


2.15 


0.27 


0.24 


0.25 


0.43 


23.1 


B+ D+(<Pir)D° 


85.17 


3.89 


3.77 




11.15 


46.4 


B+ D+{K*°K + )D° 


92.89 


5.29 


4.97 




18.66 


42.1 



D. Discussion of the results 

The values of branching fraction for the decays B + — > D~K + tt + and B + — > D+D° are 
collected in Table V. Here the 'Sj n t' error is due to uncertainties in the branching fractions 
for the decays of intermediate resonances present in the decay in question (Table I). The 
significance of the signal (Table V) is evaluated according to the formula a/— 2ln(C /C), 
where C is the likelihood calculated for the nominal fit and Cq is the respective likelihood 
function for a fit with the number of signal events fixed to zero. 

The final results for the branching fractions for the channels studied, according to Table V, 
are as follows: 

For D s — > </>(—> K + K~)n decay modes: 

B(B + -> D;K + n + ) = (1.77+° f 2 (stat) ± 0.16(sysf) ± 0.23(B int )) x 1(T 4 (5) 
B(B + -> D+D°) = (8.52t°;!(stat)) x 1CT 3 (6) 
and for D s -> i?*(892)°(^ K~ir+)K decay modes: 
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B{B + -> D-K+U+) = (2.15±g;|I(stat) ± 0.25(sys£) ± 0.43(B int )) x 10~ 4 (7) 
£(£+ -> D+D°) = (9.29t° H(stat)) x 10" 3 (8) 

The results presented here are compatible with the values reported by the BaBar collab- 
oration [9]. The branching fraction for the B + — > DfD° decay is consistent with the world 
average value given by Particle Data Group (PDG): B PDG (B+ -> D+D°) = (1.09 ± 0.27%) 
[6]. 

The same study is also performed for a large statistics MC samples simulating the fol- 
lowing processes: e + e~ — > B + B~, B°B° and qq(q = u,d,s,c). These samples are sub- 
jected to the same analysis procedure as data; the signal from the decay B + — > D^D° 
is fitted as described in subsection IV. A. The resulting branching fractions for the de- 
cays B + -> D+(-^ 07r+) J D° and B + -> Df(-^ K*°K + )D° are (9.28lg;|) x 10~ 3 and 
(8. 92I0J2) x 10~ 3 , respectively. They are compatible with the values assumed in the MC 
generator: B(B + — > DfD°) = 9.06 x 10~ 3 (Differences within statistical errors of the ex- 
tracted MC values are used as conservative estimates of systematic uncertainties for the 
reconstruction efficiencies, as described above.) 

E. Studies of two-body subsystems D S K , D s ir and Kir 

The study of two-body subsystems in the D s Kir final states is driven by two facts. The 
first one is the lack of resonances in the D + ir~ invariant mass above the value of 2.55 GeV/c 2 
observed by Belle [10], where the respective three-body decay is B + — > D + 7t~tt~. Second, 
the D~7r + and K + tt + pairs cannot form 'standard' qq resonances, so any signal observed in 
these subsystems would indicate the presence of some exotic states such as hybrid mesons, 
tetraquarks, etc. 

Below we concentrate on the D~K + n + final state. Figures 10-12 show the Dalitz plots 
of all possible combinations of invariant-mass squared: M 2 (D S K), M 2 (D s 7i) and M 2 (Kir). 
All distributions are shown both for the signal region and the AE sidebands. The structure 
visible in the invariant mass of the D s n subsystem in the mass region around 3 GeV/c 2 
corresponds to the decays B + — > (cc)K + . Preliminary studies of these Dalitz plots do not 
confirm contributions of any exotic states in the two-body subsystems. However, further 
studies might reveal some enhancements such as that observed by BaBar in the m(D~ K + ) 
spectrum [9]. Detailed analyses of the two-body subsystems are therefore still in progress. 

V. ACKNOWLEDGMENTS 

We thank the KEKB group for the excellent operation of the accelerator, the KEK 
cryogenics group for the efficient operation of the solenoid, and the KEK computer group 
and the National Institute of Informatics for valuable computing and Super-SINET network 
support. We acknowledge support from the Ministry of Education, Culture, Sports, Science, 
and Technology of Japan and the Japan Society for the Promotion of Science; the Australian 
Research Council and the Australian Department of Education, Science and Training; the 
National Science Foundation of China and the Knowledge Innovation Program of the Chinese 
Academy of Sciences under contract No. 10575109 and IHEP-U-503; the Department of 
Science and Technology of India; the BK21 program of the Ministry of Education of Korea, 



12 



the CHEP SRC program and Basic Research program (grant No. R01-2005-000-10089-0) of 
the Korea Science and Engineering Foundation, and the Pure Basic Research Group program 
of the Korea Research Foundation; the Polish State Committee for Scientific Research; 
the Ministry of Education and Science of the Russian Federation and the Russian Federal 
Agency for Atomic Energy; the Slovenian Research Agency; the Swiss National Science 
Foundation; the National Science Council and the Ministry of Education of Taiwan; and the 
U.S. Department of Energy. 



[1] Throughout this paper, the inclusion of the charge conjugate mode decay is implied unless 
otherwise stated. 

[2] S. Kurokawa and E. Kikutani, Nucl. Instr. and. Meth. A 499, 1 (2003), and other papers 

included in this volume. 
[3] A. Abashian et al. (Belle Collab.), Nucl. Instr. and Meth. A 479, 117 (2002). 
[4] Z. Natkaniec et al. (Belle SVD2 Group), Nucl. Instr. and Meth. A 560, 1 (2006). 
[5] G. C. Fox and S. Wolfram, Phys. Rev. Lett. 41, 1581 (1978). 
[6] W.-M. Yao et al. (Particle Data Group), J. Phys. G 33, 1 (2006). 
[7] F. James, CERN Program Library Long Writeup D506. 
[8] H. Albrecht et al. (ARGUS Collab.), Phys. Lett. B 241, 278 (1990). 
[9] B. Aubert et al. (BaBar Collab.), arXiv:0707.1043vl [hep-ex]. 
[10] A. Satpathy et al. (Belle Collab.), Phys. Lett. B 553 159 (2003). 



13 



dN 

dm((f>ir-) ■ (0.002 GeV/c 2 ) 



; 1 


ll 


[ i 

L ^ 

»i+^i-#+ t t'HfM*++ t+ ' 


1 

V. , 



m(07r-)[GeV/^ 



FIG. 3: Invariant mass distribution for (pTr~ pairs (points) together with the curve representing 
the fit result in the range of the D~ meson mass (red dashed line) for B + — > D~K + ir + decay. 
The signal is parameterized by two Gaussians, while the background is described by a second-order 
polynomial. The spectrum corresponds to the (AE, M^ c ) signal box. 



dN 



dm(K*{892)°K-) ■ (0.002 GeV/c 2 ) 
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FIG. 4: Invariant mass distribution for K*(892)°K- pairs (points) together with the curve repre- 
senting the fit result in the range of the Dj meson mass (red dashed line) for B + — > D~ K + tt + 
decay. The signal is parameterized by two Gaussians, while the background is described by a second- 
order polynomial. The spectrum corresponds to the (AE, M^) signal box. 
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dN 



dm{K+TT-) ■ (0.002 GeV/c 2 ) 




m{K + 7r-) [GeV/c 2 ] 



FIG. 5: Invariant mass distribution for K + n~ pairs (points) together with the curve representing 
the fit result in the range of the D° meson for the B + — > DfD° decay. The signal is parameterized 
by a single Gaussian, while the background is described by a first-order polynomial. The spectrum 
corresponds to the (AE, M^, c ) signal box. 
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FIG. 6: Distributions of the AE (top) and Mb c variables (bottom) for the decay B + — * D~ K + tt + , 
Dg — » (f)iT~ . Points correspond to data and the dashed (red) line represents results of the fit 
described in the text. 



16 



dN 

d(AE) • (0.008 GeV) 




-0.2 O.J5 -O.J -O.OS O 0.05 O.J O.I5 0.2 



A£[GeV] 



dN 

dM hc ■ (0.002 GeV/c 2 ) 




5.2 5.21 5.22 5.23 5.24 5.25 5.26 5.27 5.28 5.29 

M bc [GeV/c 2 ] 

FIG. 7: Distributions of the AE (top) and Mb c variables (bottom) for the decay B + — * D~ K + n + , 
D~ — * K*(892)°K~ . Points correspond to data and the dashed (red) line represents results of the 
fit described in the text. 
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FIG. 8: Distributions of the AE (top) and M\y C variables (bottom) for the decay B + — > DfD°, 
Df — » 4>tt + . Points correspond to data and the dashed (red) line represents results of the fit 
described in the text. 
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FIG. 9: Distributions of the AE (top) and M\y C variables (bottom) for the decay B + — > DfD°, 
Dg — > K*(892)°K + . Points correspond to data and the dashed (red) line represents results of the 
fit described in the text. 
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FIG. 10: The Dalitz plots of Mp sK vs Mp s7T for the signal box (upper plot) and AE sidebands 
(lower plot). 
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FIG. 11: The Dalitz plots of vs Mj^ for the signal box (upper plot) and AE sidebands (lower 

plot). 
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FIG. 12: The Dalitz plots of M| )< K vs for the signal box (upper plot) and AE sidebands 

(lower plot). 
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